% PhaseSpaceDensity: Computes the phase-space density
% N = number of atoms
% T = temperature (K)
% omega = trap mean angular frequency (rad/s)
function D = PhaseSpaceDensity(N, T, omega)
    cst = Constants;
    D = N * (cst.hbar * omega ./ (cst.kB * T)).^3;
end